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Supernovae connected with long-duration gamma-ray bursts (GRBs)^ are hyper- 
energetic explosions resulting from the collapse of very massive stars (~ 4OM0, where 
Mq is the mass of the Sun) stripped of their outer hydrogen and helium envelopes. 
A very massive progenitor, collapsing to a black hole, was thought to be a requirement 
for the launch of a GRB.*^ Here we report the results of modelling the spectra and 
light curve of SN2006aj,^ which demonstrate that the supernova had a much smaller 
explosion energy and ejected much less mass than the other GRB-super novae, sug- 
gesting that it was produced by a star whose initial mass was only ^ 2OM0. A star of 
this mass is expected to form a neutron star rather than a black hole when its core 
collapses. The smaller explosion energy of SN 2006aj is matched by the weakness 
and softness^'^ of GRB 060218 (an X-ray flash), and the weakness of the radio flux of 
the supernova. Our results indicate that the supernova-GRB connection extends to 
a much broader range of stellar masses than previously thought, possibly involving 
difl"erent physical mechanisms: a 'collapsar'® for the more massive stars collapsing to 
a black hole, and magnetic activity of the nascent neutron star^^ for the less massive 
stars. 
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Like all other GRB-supernovae, SN 2006aj is of type Ic.^ Its spectra resemble those 
of the dim, broad-lined, non-GRB supernova SN 2002ap.^^'^^ However, SN2006aj shows 

surprisingly weak oxygen lines for a type Ic supernova. For a comparison of the spectrum 
of SN 2006aj to those of SN 2002ap and of the GRB-supernova 1998bw, see Supplementary 
Information. 

To reproduce the spectrum of SN 2006aj^ we started from the model that was used 
for SN 2002ap,^^ but to improve the spectral fits we reduced the masses of both oxygen 
and calcium significantly, and decreased the ejected mass Mej and the kinetic energy 
accordingly. The series of synthetic spectra is shown in Figure 1. 

A lack of oxygen fines in the spectrum suggests a small Mej , but it does not necessarily 
mean absence of oxygen in the ejecta. Our model contains ~ I.3M0 of oxygen. Oxygen is 
therefore still the dominant element, but its abundance relative to other (heavier) elements 
is much lower than in SN 2002ap or in the other GRB-supernovae. Modelfing also indicates 
that oxygen is confined to high velocities (Fig. 1). A shell of oxygen comprising ~ O.IMq 
and expanding at velocities between 20,000 and 30,000 km s^^ is detected, which may be the 
result of the episode of interaction that was responsible for the early ultraviolet brightening. 

The spectroscopic results are confirmed by models of the light curve. A synthetic fight 
curve computed using the one-dimensional density and chemical abundance structure ob- 
tained from the spectral analysis reproduces the optical-infrared bolometric light curve of 
SN 2006aj (Figure 2). For SN 2006aj we derive MejRi 2Mq and 2 x 10^^ erg. These values 
are much smaller than those of the other GRB-supernova, which typically have Mejfs lOM© 
and £Jk~ 3 x 10^^ erg.^"''' The smaUer £'k and Mej involved for SN 2006aj explain why 
the light curve evolves more rapidly than that of SN 2002ap: the timescale of the light 
curve depends in fact roughly on M^^ / E^^}^ The supernova ejecta contain 0.21Mq of ^^Ni, 
which is responsible for the supernova luminosity. About O.O2M0 of this is located above 
20,000 km s~^ and causes the fast rise of the light curve. The presence of ^^Ni at high veloci- 
ties is unlikely to be the result of a spherically symmetric explosion. In a realistic aspherical 
explosion, high- velocity ^^Ni may be copiously produced near the direction of the GRB jets.^^ 

Observations in the nebular phase, when the forbidden [O l] 6300 and 6363 A lines should 
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be strong in emission, will be needed to determine more accurately the value of Mej. Such 
observations, to be performed starting August 2006, will also be useful to study any possible 
asymmetry and the orientation of the supernova with respect to the line of sight to the 
Earth, and thus to link the supernova with the GRB.^^'^^ 

The properties of both the supernova (small energy, small ejected mass, low oxygen 
content) and those of the GRB (unusually soft and long) seem to suggest that GRB 060218- 
SN 2006aj was not the same type of event as the other GRB-supernovae known thus far. The 
radio properties of SN 2006aj were also intermediate between those of the GRB-supernovae 
and of SN 2002ap.ii 

One possibility is that the initial mass of the progenitor star was significantly smaller 
than in the other GRB-supernovae, and that the collapse/explosion generated less energy. 
A star of zero-age main-sequence mass ~ 20 — 25Mq would be at the boundary between 
collapse to a black hole or to a neutron star.^^ If the star collapsed only to a neutron 
star, more core material would be available to synthesize ^^Ni. For example, a star with 
~ 20Mo initially would have a carbon-oxygen core of ~ 3.3Mq}^ If core collapse left behind 
a neutron star of ~ IAMq, ~ 1.3M0 of oxygen and ~ O.6M0 of heavier elements (including 
both intermediate-mass elements such as Si and Fe-group elements) could be ejected in 
the supernova, consistent with our results. Such a collapse is thought to give rise to an 
explosion of ^'k^ 10^^ erg, but there are indications of a spread in both and the 
mass of ^^Ni synthesized. ^° Additionally, magnetar-type activity may have been present, 
increasing the explosion energy. Magnetic activity may also have caused the very long 
duration of the 7-ray emission^^ and the mixing-out of ^^Ni required by the rapid rise of the 
hght curve. It is also possible that in this weaker explosion the fraction of energy channelled 
to relativistic ejecta was smaller than in the classical GRB-supernovae, giving rise to an 
X-ray flash (XRF).^! 

Another case of a supernova associated with an XRF has been reported. The puta- 
tive supernova, although poorly observed, was also best consistent with the properties of 
SN 2002ap.^^ This may suggest that XRFs are associated with less massive progenitor stars 
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than those of canonical GRBs, and that the two groups may be differentiated by the forma- 
tion of a magnetar^^ or a black hole, respectively. The properties of both the GRB and the 
supernova may scale with the mass of the progenitor. Still, the progenitor of SN 2006aj 
had been thoroughly stripped of its H and He envelopes. This is a general property of all 
GRB-supernovae known so far, and possibly a requirement for the emission of a high energy 
transient, which may be more easily achieved in a binary system. ^^'^^'^^ 

If the star was initially more massive (> 25M0), and it collapsed directly to a black hole 
as in the more powerful GRB-supernovae events, a number of questions arise. Why was 
the energy of the explosion so small? Where did the large core mass end up? Continuing 
accretion onto the black hole could explain the missing mass. This may occur if the angular 
momentum of the core was smaller than in the more energetic cases. Other more exotic 
scenarios, such as merger models, may also work. 

A case of a progenitor mass just exceeding the black hole limit may be that of SN 2002ap. 
This SN may not have produced a magnetar and an XRF, because it did not collapse to a 
neutron star but rather to a black hole,^^ yet at the same time the energies involved in the 
collapse may have been too small to give rise to a GRB. 

In our scenario, some soft 7-ray repeaters energized by a magnetar^^'^^ may be remnants 
of GRB 060218-like events. Magnetars could thus generate a GRB at different times. As 
they are born, when they have a very short spin period (~ 1 ms), an XRF (or a soft GRB) 
may be produced as in SN2006aj-GRB 060218. Later (after more than 1,000 years), when 
their spin rate is much lower, they could produce short-hard GRBs by a giant flare.^^ Finally, 
if the progenitor star had a massive companion in a close binary system, as may be required 
for the outer envelope to be stripped and a long-duration GRB or XRF to be produced, 
the system may evolve to a close double-neutron-star system. When the two neutron stars 
finally merge, a short-hard GRB may again be produced. 
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Figure 1. Spectra of SN2006aj and synthetic fits. The observed spectra of SN2006aj 
(coloured traces) are cahbrated in the V band, but elsewhere they my be distorted,^ hence the 
poorer agreement in some of the red parts. Also, the blue part is not reliable shortward of ~ 4200 A. 
The synthetic spectra (black traces) were computed using our Monte Carlo spectrum synthesis 
code.^'^ Because of the spectroscopic and photometric similarity to SN 2002ap,^^ we used a similar 
model of the explosion, but in order to improve the match we reduced the masses of both oxygen 
and calcium significantly, and decreased Mgj and accordingly. Our model has Mej~ 2Mq and 
2 X 10^^ erg. The strongest features in the spectra are due to lines of Fe ii, Ti ii, and in the 
later phases Ca ii (< 4500 A), Fe iii and Fe ii (near 5000 A), Si ii (near 6000 A), O i (near 7500 A), 
and Ca II (near 8000 A). The O i and Ca ii lines become stronger at more advanced epochs, and 
are conspicuous because they form at a roughly constant wavelength, corresponding to a velocity 
(~ 25, 000 km s~^) higher than that of other lines. This indicates the presence of a shell of material, 
dominated by oxygen, at velocities between about 20,000 and 25,000 km s~^. This high- velocity 
material may result from the piling up of circumstellar material on the expanding ejecta. We mod- 
elled the spectrum by adding a small amount of mass (~ O.IOM©) at 20, 000 <v< 30, 000 km s^^. 
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Figure 2. The light curve of SN 2006aj. The bolometric Hght curve of SN 2006aj (circles) 
is compared with the model light curve (solid line), and with the bolometric light curve of SN 
2002ap (squares). A supernova light curve is powered by 7-rays released in radioactive decays 
of freshly synthesized unstable ^^Ni to ^^Co and hence to stable ^^Fe. The 7-rays deposit in the 
dense cjccta, giving rise to a flux of optical photons. The light curve rises at first as the diffusion 
time of photons decreases as the ejecta expand. A maximum is reached when the escaping photon 
luminosity approximately equals the deposited energy. -"^^ The light curve then declines as the 
density becomes low enough to allow significant 7-ray escape. The more massive the supernova 
ejecta and the smaller their kinetic energy, the more difficult it is for photons to escape, which 
means that the light curve reaches maximum later and has a broader peak. The bolometric light 
curves were constructed by integrating the optical and near-infrared fluxes (for SN 2006aj, optical 
photometry obtained with European Southern Observatory's (ESO) Very Large Telescope (VLT) 
and ncar-infrarcd photometry reported in the Gamma-Ray Burst Coordinates Network (GCN) 
were used), after correcting for the host-galaxy distance/redshift and the reddening toward the 
supernova - for SN 2006aj, 143 Mpc, z = 0.0335, and E{B - V) = 0.126 mag'^]. The model light 
curve is synthesized from the one-dimensional density and chemical abundance structure of the 
best-fitting spectral models. It corresponds to ~ 2Mq ejecta expanding with a kinetic energy of 
~ 2 X 10^-*^ erg, having in total 0.2Mq of ^^Ni. The small amount of mass and energy added by 
the inclusion of the outer oxygen shell (see Fig. 1) have a very limited impact on the light curve 
since the mass is located at low densitv and has low ootical deoth. The exnlosion of SN 2006ai is 



